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ABSTRACT 

We present an upper limit for the relative escape fraction (/JJ^^) ionizing 
radiation at z~3.3 using a sample of 11 Lyman Break Galaxies (LBGs) with deep 
imaging in the U band obtained with the Large Binocular Camera, mounted on 
the prime focus of the Large Binocular Telescope. We selected 11 LBGs with 
secure redshift in the range 3.27 < z < 3.35, from 3 independent fields. We 
stacked the images of our sources in the R and U band, which correspond to 
an effective rest-frame wavelength of 1500A and 900A respectively, obtaining a 
limit in the U band image of >30.7(AB)mag at la. We derive a la upper limit of 
fefc ~ 5%, which is one of the lowest values found in the literature so far at z~3.3. 
Assuming that the upper limit for the escape fraction that we derived from our 
sample holds for all galaxies at this redshift, the hydrogen ionization rate that we 
obtain (r_i2 < 0.3s~^) is not enough to keep the IGM ionized and a substantial 
contribution to the UV background by faint AGNs is required. Since our sample 
is clearly still limited in size, larger z~3 LBG samples, at similar or even greater 
depths are necessary to confirm these results on a more firm statistical basis. 

Subject headings: galaxies: distances and redshifts — Galaxies: evolution — 
Galaxies: high-redshift — intergalactic medium — diffuse radiation 

1. Introduction 

The evolution of the intergalactic medium (IGM) depends primarily on its ionization 
state, which at high redshifts is a function of the ionizing UV background (UVB) produced 
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by galaxies and AGNs. In the early universe, star forming (SF) galaxies are ubiquitous 
and in principle can provide the necessary flux to keep the IGM ionized at least up to 
z ~ 6, when the re-ionization appears to be completed (jPan et al. 2006 ). However, the 
poor knowledge of the average escape fraction of ionizing Lyman continuum photons from 
the interstellar medium of each galaxy introduces large uncertainties to the exact level of 
contribution in ionizing the IGM. The escape fraction of AGN is not strongly variable with 
redshift but the apparent number density of bright QSOs and AGNs is rapidly decreasing 
at 2 > 3, thus it is assumed that the contribution to the ionizing f lux of the SF gala xies 



should become dominant at 2; > 3 (e.g. iFaucher-Giguere et al.ll2008l : ICowie et al.ll2009l ). In 



this respect, several attempts have been made to derive the escape fraction of UV ionizing 
photons both at low redshifts, from space, and at high redshifts [z = 2—4) from ground based 
observations. In most cases only upper limits on the escape fraction were obtained, giving 
little evidence to support a scenario where enough ionizing photons escape from galaxies. 
Actually, early measurements at low and intermediate redshifts, with fesc <l-3%, seem to 
indicate that galaxies are not the major contributo rs to the ionizing UV background at z < 2 



leaving AGNs as t he main ionizing population (e.glGia 



Siana et al.ll2007t iCowie et all 12010 : iSiana et al.ll2010l : 



longo et al. 



Bridge et al. 



997l:lMalkan et al.ll2003l : 



20101 ). 



At high redshifts direct observations of the Lyman continuum flux for galaxies at z > 4 
are difficult because of the sharp increase in the neutral hydrogen absorption by the IGM. 
For ground-based surveys, the search has thus been focused at 2; ~ 3 — 4 which is a com- 
promise between lower detection efficiency in the UV and higher transparency of the IGM 
in comparison to higher z. In this respect, there are two main strategies to face the issue 
of escaping ionizing photons from galaxies. One way is to derive direct constraints using 
very deep UV imaging and/or spectroscopic observations in order to ascertain if even a 
small fraction of LyC escapes from typical star forming galaxies. Such surveys are time 
consuming, requiring 20-30 hours of integration for reasonable UV efficiency, and for this 
reason inv olve small samples rn ade up of some tens of galaxies. T he works by 



Steidel et al. 



( 1200 ll ) and lShapley et al.l (|2006[ ) are typical in this respect. Indeed, ISteidel et al.l (l200l[ ) first 
claimed appreciable ionizing flux corresponding to a relative (with respect to the observed 
I5OOA flux) escape fraction of the order of > 50% from a composite spectrum of the bluest 
29 galaxies of their large LBG sample at z ~ 3.4. If such an escape fraction were typical 
of the LBG population, this would imply an ionizing UVB at least 3 times higher than the 
one derived from t he analysis and mo d eling of the Lyrnan a a bsorption by the IGM at the 
same redshift (e.g lBolton et al.l 120051 ). iGiallongo et al.l (120021 ) obtained the flrst high S/N 
long slit spectra of two LBGs at z = 3 and z = 3.3, selected from Steidel's sample, where 
no signiflcant Lyman continuum emission has been detected. This implies a la upper limit 
to the relative escape fraction of 15%. Using individual spectra of 14 LBG galaxies at z ^ 3 
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Shapley et al.l ( 120061 ) found significant emission only in two galaxies, implying a detected 
average relative escape fraction of 14% for their sample at a 3 a confidence level, which is 
~4.5 times lower than the value derived by ISteidel et al.l (l200l[ ). The fact that the estimate 
of the escape fraction value obtained in this way has been progressively reducing from 70% 
to 15%, reveals the difficulties and the uncertainties involved in such measurements. 

Another way to measure escaping ionizing flux is to use shallower UV imaging but of a 
very large sample of galaxies. In this case, instead of measuring the average/typical escape 
fraction of the population, the effort is concentrated in finding possible small fractions of 
special types of galaxies which could show a large percentage (>50%) of escaping ionizing 
flux. The main limitations of this method are the absence of extensive spectroscopic confir- 
mation of the redshifts in the sample or, whenever spectral information is available, the large 
redshift range of the sources, which requires extended simulations to correctly account for 
the variances of IGM absorption, and the lower S/N leve l of the narrow band (NB) images. 
Typical in this direction is the work by llwata et al.l (120091 ). who have used a UV narrow band 
filter technique and found ionizing radiation, at a 2a level, in 7 out of 73 LBGs at 2; > 3 
LBGs and 10 out of 125 Lya emitters at 2 > 3 out of a sample of 198 z > 3 galaxies in the 
SSA22 field. The relative escape fraction derived for these 7 detected LBGs reaches 80%. 
It is interestirig to ri ote that they report null emission for the object SSA22a-D3, for which 
Shapley et al.l (120061 ) claimed an ionizing fiux dete ction. Moreover in several cases (including 
the second object detected by lShapley et al.ll2006l ). an offset is present between the position 
of the object in the UV narrow band filter and in th e R ban d detection image, suggesting 
possible contamination by interlopers. IVanzella et al.l (l2010bl ) evaluated the probable cause 
of such offsets, showing that there is a > 50% probability for contamination, by lower z 
objects, in at l e ast 30% of the galaxies showing clear Lyman continuum detection. Recently, 
Vanzella et al.l (l2010al ) using ultra-deep VIMOS intermediate U band and deep FORSl nar- 
row band imaging of 102 galaxies in the GOODS-South at z ~3.7, obtained an f^sc upper 
limit of 5%-20% at a 3a level, depending on the assumed extinction curve and E(B-V) value, 
as well as the luminosity and redshift of the sources adopted in the analysis. 

These two approaches should be considered as complementary, since they give us infor- 
mation on different aspects of the problem: the typical escape fraction from the majority 
of galaxies, the former, and the existence of rare galaxies with higher escape fractions, the 
latter. Both methods are necessary to measure with greater accuracy the global amount of 
ionizing flux in the universe. In the present paper we use deep U and R ba nd imaging ob- 
tained by the UV optimized Large Binocular Camera ( Giallongo et ahlbood ) in two Steidel 
flelds and in a region of the COSMOS fleld to derive stringent limits on the ionizing escape 
fraction of 11 LBGs at z ^ 3.3. The uniqueness of our data is represented by the simultane- 
ous availability of very deep UV images, obtained with one of the most efficient ground based 
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large field UV imager at an 8m class telescope, in fields where spectroscopically confirmed 
galaxies in a very narrow redshift range around z ~ 3.3, are already available. This sample, 
although it is limited in size, gives us the possibility to focus on a very small redshift range, 
reducing appreciably the effect of the variation with redshift of the IGM opacity and allowing 
a more direct measure of the relative ionizing escape fraction from LBGs at 2; ~ 3.3. 



Data 



The targets were selected from three fields where deep UV images we re obtained wit h 
the Large Binocular Camera (LBC) at the Large Binocular Telescope (LBT. lHill et al.ll2010l ). 
The LBC has two channels, one optimized for the UV and blue bands and one optimized 
for imaging in the red bands, that was used to obtain the R band images. The detectors 
have four 4Kx2K chips with a pixel scale of 0.23 arcsec/ pixel, pro viding an unvignetted 
field of view of about 23 x 23arcmin'^ . The standard LBC pipeline ( iGiallongo et al.ll2008[ ) 
has been used for data reduction. In particular, after creating the bias-subtracted images, 
we have applied standard fiat-fielding using s ky fiats obtained durin g twilight. The sky 
background was subtracted using SExtractor (IBertin &: Arnoutd Il996[ ) with a background 
mesh of 64x64 pixels and adopting a median filter of 3x3. The astrometric solution was 



com puted using t he software AstromC developed by iRadovich et al.l (|2004[ ) and the USNO- 
A2 ( iMonet Ill998l ) as reference catalog. After applying the astrometric correction, we have 
coadded the various images to create deep mosaics for each field. We did not include a color 
term in the photometric calibration and the accuracy of the zeropoint is 0.03-0.05 for the U 
filter. P art of these irnages h ave already been used to derive the number counts in the U 
band by iGrazian et al.l ( l2009l ) providing good agreement with previous results. 



For each scientific image the LBC pipeline computes a corresponding RMS map directly 
from the raw science frames. We obtained the absolute RMS maps for each pointing by com- 
puting the RMS in each individual image (using Poisson statistics and the instrumental gain 
and self-consistently propagating this RMS over the whole data reduction 



as cr,' 



Nadu 
gain 



process. The delivered RMS maps were found to be consistent with the absolute RMS ex- 
pected for each image a nd with the median R MS derived from 5x5 pixel boxes in random 
positions of the sky (see iGawiser et al.l (120061 ) for details). 



Two LBC fields come from Steidel's sample and include the bright quasars Q0933-I-28 



fa=09:33:31 S=+ 28:U-A2) and Q1625+26 (a=16:25:30, 5=+26:52:31) ISteidel et al.ll2003l 
Reddy et al.ll2008l ). A major advantage that led t o the selection of the se areas is the existence 
of several spectroscopic redshifts, obtained by ISteidel et al.l ( 120031 ) in the redshift range 
2 < z < 3.5. The LBC images were obtained in the period 2007-2009 in the U band with 
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total exposure times of 12h and 8h, respectively, and in the R band with total exposure 
times of 2.8h and 3.4h. The AB magnitude limits reached at Icr in the U band are 28.8 and 
28.6 and in the R band 25.5 and 25.9 at lOtx, respectively, using an aperture of 2xFWHM. 
The seeing in the R band is 1.0" for Q0933+28 and 0.8" for Q1625+26. For both fields, the 
average seeing in the U band is 1.1". 



The third LBC field is within the COSMOS field flScoville et all 120071 : iLillv et all 120071 ) 

where extensive multiband and spectroscopic dataset is continuously improving. The LBC 
observations of the COSMOS field started in February 2007 to add deep UGRIZ imaging in 
the central 1 deg^ and are currently on-going. The exposure time per pointing in the U band 
is ~6 hours, for a total of three partially overlapping LBC pointings, reaching 11 hours at 
the deepest region. Actually, the whole area covered by deep UGRIZ data in COSMOS is 
~700 sq. arcmin. and reaches a la AB magnitude limit of U=28.7, with a typical seeing of 
1.0". The R band image has an average seeing of 0.9" and it reaches an AB magnitude limit 
of R=25.9 at lOcr, for two pointings. The third pointing was obtained during commissioning 
and it has a total exposure time of 0.5 hours, an average seeing of 0.9" and an AB magnitude 
limit of R=23.8. The datasets in all three fields are composed by different mosaics, obtained 
in different years and in a variety of seeing conditions and exposure times. Although we have 
average seeing values estimated in each band, these values do not remain stable throughout 
the whole final mosaic. Actually, the seeing in the R band ranges from 0.8" to 1.0" and in 
the U band from 1.0" to 1.2", especially in the COSMOS field where we cover the largest 
area and we have 3 different pointings with different final exposure times. 



3. Analysis 



It is important to note that there are several definitions for the ionizing escape frac- 
tion which depend on the non-ionizing flux considered as a reference. The absolute escape 
fraction is the ratio between the escaping Lyman coii tinuum flux and the one intrinsically 
produced by stars in the galaxy (ILeitherer et al.lll995[ ). In practice this value is difficult to 
determine because we should know the intrinsic Lyman continuum flux from an accurate fit 
to the overall spectral energy distribution of the galaxy. For this reason, observationally, 
the ionizing escape fraction is usually related to the rest frame 1500A flux and denoted as 
fesc- This reference flux is however attenuated by dust and a correction is needed if we want 
to derive the intrinsic Lyman continuum. If the amount of dust attenuation is not known 
in the analyzed galaxy sample, a relative escape fraction f^^j: can be introduced as the frac- 
tion of escaping Lyraan c ontinuum photons divided by the fraction of escaping photons at 



1500A ( Steidel et al.ll2001 ). The relation between the two quantities is fesc = fesc ^ 10" 



0.4Ai 



- 6 - 



where Aisoo is the dust absorption at 1500A ( Inoue et aL 2006 : Vanzella et al.ll2010a ). Since 
LBGs at z=3 show an average ^1500 =0.6 flVanzella et al.ll2010ah then /J,^^ should be typi- 
cally 2 times larger than f^sc- Finally, to derive the Lyman continuum ionizing fraction from 
the observed fluxes at 900A and 1500A , we need to estimate the contribution to the ab- 
sorption by the intervening IGM and the average intrinsic UV spectral shape of the galaxy 
populations. In our analysis we estimated the relative escape fraction using the following 
equation: 



f 



rel 
esc 



(-^1500/ Lgoo)int ,_r_IGM\ 



(/l50o/ /900) ohs 



■exp(r< 



900 



where (Ivi5oo/-^v9oo)mt is the average intrinsic ratio of non-ionizing to ionizin g specific in- 
tensities as derived from spectral synthesis models (e.g. Bruzual fc CharlotlboOS ). (/1500//900) obs 



is the observed flux ratio estimated from our U-R color and exp{TgQQ^) is the inverse of the 



average IGM transmission at 2; ~ 3.3. For comparison with previous estimates of f^gj: in the 
literature, we adopt (Li5oo/i^9oo)jnt=3, but this value is model dependent and will be dis- 
cussed below. A major issue in the derivation of the relative escape fraction is the estimate 
of the average attenuation produced by the IGM along the line of sight, exp^r^QQ^'^). This 
value is very sensitive to the actual redshift distribution of the sources and to the wavelength 
bandpass of the filters used. The broader the filter, the more the luminosity of the continuum 
at A < 912A will be diluted, resulting to a less constraining value for the escape fraction. 

The average IG M transmission expf— ri m^"^) for our galaxy sample has been derived 
following the work of iProchaska et al. where they analyzed the average rest-frame 

spectra of z>3 SDSS quasars. They measured the IGM absorption from the drop of the flux 
at A <912A, rest frame, in a redshift range 3.6< z <4.3. To reproduce the IGM attenuation 
at z ~3.3, instead of using the Lyman absorption statistical distribution in column density 
and redshift, we decided to adopt the same empirical fit described by the authors. More 
specifically, since the lowest redshift they consider is 3.6, we extrapolated the IGM LyC 
absorption to th e avera ge redshift of our sample (z=3.3) using the formula (Eq.6) given by 
Prochaska et al.l (120091 ). a nd adding the con tribution of the Lyman series absorption using 
the equations described by iFan et al.l (120061 ) . The attenuation of photons that were emitted 
at an observed wavelength, Xobs, at redshift 2em=3.3, is shown in Fig. 1 (solid black line). 

The LEG UV broad band transmission covers the wavelength range from 32OO-39OOA, 
with a small tail up to 4000A (considering an average airmass for our dataset of 1.2 and 
including instrument and GCD response). This means that Lyman continuum emission by 
LBGs at A <912A starts just short-wards of the red edge of the U filter for sources with 
z ~ 3.3, leaving the LBGs non ionizing emission outside the filter at Xobs > 4000A . The 
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effective rest-frame wavelength for the Prochaska IGM attenuation, convolved with the LBC 
U band filter and considering z=3.3, is 86OA , which is reasonably close to 900A . Thus, in 
order to avoid contamination from continuum emission red-wards of the Lyman limit and to 
probe a wavelength close to 900A , without introducing excessive dilution of the observed 
emission due to IGM attenuation, we decide to use only LBGs with spectroscopic redshifts 
in the range 3.27 < z < 3.35. We consider this as the best trade-off in order to obtain tight 
constraints for the escape fraction. 

We have a total of 12 LBGs in the tight redshift range (3.27 < z < 3.35) in all of 
our fields: 4 sources in Q0933, 1 in Q1623 and 7 in COSMOS. In the COSMOS field one 
of the sources (a=10:01: 19.29, 5=+02:04:20.2) showed significant emission in the U band 
(> 3(t), which was slightly offset from the image center in the R band. We checked the high 
resolution ACS image finding that the morphology of the source appears more consistent 
with two distinct objects, which are too close to be resolved in ground based observations. 
Careful inspection of the spectrum also supports such a scenario, where the main spectral 
features clearly appear only in the UV dropout candidate. In the LBC images there is an 
offset of ~ 0.6", between the emission in the R and U band that corresponds to ~6kpc. 
For all these reasons, we decided not to inclu de this source in our sa mple, leaving us with a 



total of 11 LBGs in the 3 fields. According to lVanzella et al.l ( l2010bl ). assuming the number 



counts for U<28.5 with a seeing of 1", there is more than 80% probability that 10% of the 
sources are affected by contamination from lower redshift sources in the line of sight. Thus, 
the fact that 1 out of our 12 LBGs seems to be contaminated is within the statistical limits. 

To compute a reliable upper limit to the Lyman continuum escape fraction we created 
thumbnail images of 45" x 45" in the U and R bands around each of the selected sources. 
The size of the thumbnail is large enough to obtain a reliable estimate for the background, 
which is subtracted after masking the surrounding objects. We then sum all thumbnails in 
each band creating a stacked image, where the flux corresponds to the weighted mean value 
of the fluxes of the individual images in each band. In Fig. |2] we show the result of this 
procedure. We obtained aperture photometry around the source in the R band (/1500) where 
the source is visible, extracting the center of the stacked profile. This position is then used to 
measure the expected aperture flux in the U band image (/goo)- Fig. [2] shows no significant 
flux in the U band stacking. In order to establish the optimum apertures for the photometry 
we created a stacking of stars. The selected stars are near our candidates in the 3 fields, 
and their fiux has been normalized in order to get the same magnitude distribution of the 
LBG stacking. This is necessary since the magnitudes of the LBGs used vary significantly 
and the characteristics of the brightest sources will be more prominent in the final stack. 
The average measured stellar FWHM is 0.9" in the R and 1.1" in the U band. Using the 
stellar PSF, we have created a kernel that matches the PSF in the R band to the one in 
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Fig. 1. — The IGM transmission is calculated using the iProchaska et al.l ( 120091 ) value ex- 
trapolated at z=3.3 (black solid line) along with the total transmission curve for our U band 
(red solid line). The lower black dotted line corresponds to the Madau IGM transmission 
and it is shown for comparison only. The solid and dotted lines under the filter transmission 



are the resulting exp{- 



J-IGM-] 
'900 ) 



curves after the convolution with the filter response using the 



Prochaska and Madau IGM transmission curves respectively. The color-coded lines corre- 
spond to the upper and lower limit o f the IGM attenuation showing the uncertainty in our 
measurement as given by the work of IProchaska et al.l (120091 ) (see parag. 4 for details). The 
filter transmission and the resulting exp{—Tg^^) curves after the convolution have all been 
multiplied by a factor of 2 to facilitate the view. At A ~ 3700 the real filter transmission 
value is 0.19. 
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the U band. Once calculated, such Gaussian filter has been applied to the R band image of 
the LBG stacking, in order to match the R band morphology to the expected U band size. 
The finally adopted aperture diameters correspond to l.SxFWHM in U band, namely, 7.3 
pixels. The magnitude measured in the smoothed R band is 24.35 mag (AB). The la rms 
background error in the U band has been used as the flux upper limit and it is equivalent to 
a lower limit in magnitude of U > 30.7(^45) at la confidence level for the LBG stacking. 



4. Results and Discussion 



Using the 11 sources with spectroscopic redshift z~3.3 we obtain (/i5oo//90o)ofes >224. 
Considering the average IGM attenuation extrapolated at this redshift from the Prochaska 
value (exp(— rg(^*^)=0.27), we derive fj^^^ <0.05 (5%) at la. Compared to the values in 
the literature so far, this is one of the most constraining results for this redshift range. 
An additional advantage of our sample is that it covers 3 independent areas in the sky 
and includes spectroscopically confirmed LBGs that have been selected through different 
methods. More specifically, the LE iGs in the Q0933 f ield are U drop-out sources, the one in 



the Q1623 field is a BX candidate ( iReddy et al.ll2008[ ) and out of the 6 sources in COSMOS 



4 have been selected as both U drop-out and gzK candidates, while 2 are only U drop-out 
candidates. 

For a more detailed approach, we used the same technique to estimate the escape fraction 
also for individual sources. The resulting values are presented in Table [H The average 
magnitude of the sources is ~24.4 in the R band and the values of the upper limits of the 
relative escape fraction for individual sources range from <6% to <36%, with an average 
value of <21%. To exclude the possibility that our result is biased by the wide range in 




Fig. 2. — The stacking result for the 11 LBGs of our sample. The size of each thumbnail is 

45" X 45" 
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the R magnitude we created a subset of our sample, where we used for the stacked image 
only the sources with 23.5 < R < 24.7 (6 sources). In this case we obtain a relative escape 
fraction of <6%. 

As already stated in the introduction, there are two main uncertainties involved in the 
measurement of an accurate escape fraction, concerning the transmission of the IGM and 
the intrinsic Li 5nn/-^flnn ratio. The u ncertainty associated to the IGM is computed following 



the work of iProchaska et al.l (120091 ). Using the la limits they state for the IGM attenua- 



tion, we compute upper and lower limits for the value of exp^TgQQ^). This way, we obtain 



upper-0-^7 and exp{-Tl^Q^^)iower=0-17, from which we derive {f^sc)iower <0.04 
p*^^'^ <0.08 



900 

(4%) and {f^gc)upper <0.08 (8%) respectively. For comparison with early results in the liter- 
ature we derive also the relative escape fraction using the Madau model. As seen in Fig.l 
the exp{-Tl^Q^') 

Madau IS 0.12 and we obtain f^f^ <0.11 (11%). Si nce Madau ' s first analysis 



the statistics of IGM abso rption have changed appreciably (e.g iFan et al.l l2006l : iMeiksin 



20061: llnoue fc Iwata I 20081 ) and Madau's model represents now a lower limit to exp(— r. 



IGM\ 
900 ) 

(jinoue fc Iwata II2008I ). making the flf^ value based on this model an upper bound. In fact, 
considering all uncertainties we find that the highest upper limit for the escape fraction in 
no case exceeds 11%, that corresponds to twice our fiducial value. 

Another source of uncertainty in estimating the relative escape fraction is the assumed 
value for the parameter Li5nn/Lflno- So far, the most wide ly used value is Li5oo/L9oo=3 (e.g. 



Steidel et al.l l200ll : llnoue et al.l l2005l : IShapley et al.l 120061 ) and we have adopted it in order 



to facilitate the comparison of our results with previous works in the literature. Such value 
is derived considering typical models of UV SED for star-forming galaxies with constant 
star formation rate (SFR). However, this ratio of the intrinsic luminosities depends also 
on the age o f the galaxy and on its star formation history (SFH). Even for constant SFH, 
Inoue et al.l (120051 ) find that the intrinsic ratio depends on the time passed since the onset 
of star formation and ranges fro m 1.5 to 5.5, approaching the upper edge for older ages. 
More specificallv Meiksin I ( 120091 ) has found a /j.(1500A )/f^(912A ) that ranges from 2.7 
considering IMyr starbursts to 6.2-6.7 for 100-300Myr that are the typical a ges of LBGs. If 



we cons ider an intrinsic luminosity ratio of 7, that is the value recently used by lVanzella et al. 



(l2010al ). our estimated escape fraction would increase by a factor of 2.3 {fill <11.5%). 



Nevertheless, assessing the role of LBGs in producing the ionizing U VB does not depen d 
on the assumed value for the parameter L1500 / L900 • Following the work by llnoue et al.l (120061 ) , 
where they define the Lyman continuum to UV escape flux density ratio, the luminosity 
density pgoo.esc corresponds to: 



P900,e 



-^900 r 

P15007 /( 



rel 
esc 



Pl500( 



^1500 



900 



1500 



lofes6Xp(TgQg ) 



(2) 
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After converting; the luminosity density uncorrected for dust extinction derived by lReddy fc Steidel 
fl2009h at 1700A to 1500A , we obtain pisoo = (3.61 ± 0.24) x W^ergs~^Hz-^Mpc-^ , 



which corresponds to a pgoo,. 



5.98 X 10 er^fs Hz ^Mpc We note that this lumi- 



nosity density has been obtained using the slope of the luminosity function (LF) given by 
Reddy fc Steidel I (120091 ) . a = —1.73 ± 0.13, and integrating this LF down to their faint 
magnitude limit (Ma_b(1500)=-17.6 corresponding to R~28), assuming that the fesc does 
not evolve with luminosity. If we consider sources with R<25.5, that is the magnitude limit 



-20.2), we obtain pi5oo = 1.27 x 10^^ erg s'^Hz-^Mpc 



^ , fr om 



of our sample (M^s(1500) = 
which we derive p9oo,esc = 2.1 x 10'^^ergs^^Hz^^Mpc^^\ At 2; ~ 3 Bolton et al.l ( 120051 ). 
using high resolution simulations to model the Lya forest opacity, estimate the hydrogen 



ionization rate of the IGM as r_ 



12 



(ionization rate in 10 s 



~^atom ^] 



which corresponds to a luminosi t y den sity of pgoo = 16.5 x lO'^^ergs ^Hz ^Mpc ^. The 
values presented by iBolton et al.l ( 120051 ) are in good agreement with more recent works (e.g. 



DairAglio et al. 20091 ). Our estimated contribution, derived from Mab(1500) <-20 sources, 
is significantly lower than the theoretical value. Even considering the contribution of fainter 
galaxies, the estimated emissivity does not increase enough in order to supply the required 
ionizing radiation. Due to our tight constraint on the escape fraction, we derive an up- 
per limit to the hydrogen ionization r ate r_i2 < 0.3s^^, using as an average spectral index 
auv=^-8 jFaucher-Giguere et al. 2008 ) and Mab(1500) <-17.6. This is lower than the values 
reported by Bolton et al. ( 20051 ) and only including their estimated AGN contribution at z~3 
(r_i2 = 0.4s~^), the ionization rate is barely consistent with the theoretical expectations 



from the Lya forest analysis. 

Here we should discuss the possibility of a bimodal distribution of the LyC escape 
fraction, where a small fraction of the galaxy population at a given redshift can be the 
major contributor to the cosmological ionizing flux. At present, our sample of 11 LBGs is 
still statistically modest and cannot exclude in a conclusive way the existence of a limited 
sub-group of galaxies with significant escape fraction. Moreover, the search for such rare 
strong LyC emitte rs at high redshifts is not as straightforward as at lower redshifts. Indeed, 
Cowie et al.l (120091 ) found only 1 galaxy out o f 600 present i ng a s ignificant fraction of escaping 
ionizing flux at z ~ 1. At higher redshifts, llwata et al.l (120091 ) claims that up to ~10% of 
galaxies at 2; ~ 3.1 show escape fractions at the level of ~ 80%. However, an offset (1" — 2") 
between the image in the R band, sampling the non- ionizing flux at ~1500A and the UV 
images, sampling the LyC flux, has been noted by the same authors. Since the escape 
fraction derived by the UV-R color in the detected LyC emitters is at a the level of ~ 80%, 
this would imply an unrealistic fraction, > 100%, from the limited spatial region where 
the claimed L yC emission is observe d. Through statistical analysis, this number has been 
questioned by lVanzella et al.l (l2010bl ). claiming that at least 30% of those galaxies have more 
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than 50 % probability of being low redshift interlopers. On the other hand, IVanzella et al. 



( l2010af ) have found only 1 galaxy out of 102 LBGs at 3.4 < z < 4.5 that shows strong 
LyC emission. In any case, assuming a population of strong LyC emitters were present at 
z=3, this could be at most of the order of 1-7% for the overall population. Adopting this 
contamination fraction as a reference value, the contribution of this subsample to the UV 
emissivity translates to an increase in the ionization rate of the order of r_i2=0.1-0.3s~^, to 
be added to our original upper limit (r_i2 < 0.3s~^). Even in this case, the ionization rate is 
still not sufficient to keep the IGM ionized, making necessary an appreciable contribution by 
some other population, possibly AGN or a large number of very faint star forming galaxies 
with increasingly larger escape fractions. 

Such a sce nario, where the A GN contribution is still important at 2; > 3 is also supported 
by the work of ISiana et al. J2008h . who derive an ionization rate of r_i2 = 0.48s ^ for QSOs 



aX z 3.2, suggesting that AGN and LBGs are emitting comparable fracti ons of ionizing 



flux to the intergalactic medium. This is in agreement with recent results by lGlikman et al. 



(120111 ). who find that the AGN account for at least half the ionizing radiation needed at 
these redshifts (60±40%). Moreover, our observed limit for the /J^^, is closer to the values 
reported by large surveys at lower redshifts {1 < z <2) suggesting a milder evolution with 
redshift. Deeper UV and narrow band observations of larger statistical samples of LBGs at 
z~3 will allow us to clarify this issue. 
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Graham, Arizona. The LBT is an international collaboration among institutions in the 
United States, Italy and Germany. LBT Corporation partners are: The University of Ari- 
zona on behalf of the Arizona university system; Istituto Nazionale di Astrofisica, Italy; 
LBT Beteiligungsgesellschaft, Germany, representing the Max-Planck Society, the Astro- 
physical Institute Potsdam, and Heidelberg University; The Ohio State University, and The 
Research Corporation, on behalf of The University of Notre Dame, University of Minnesota 
and University of Virginia. 
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Table 1: Summary of the f^f^ values for the individual LBGs. 



ID 


RA 


DEC 




z 


Rmag 


Umag 


frel 
J esc 












±0.07 


1(7 (u.l) 


la (u.l) 


3400 


143.35424 


+28.80694 


3. 


,27 


24.88 


29.75 


0.203 


12646 


143.32868 


+28.71913 


3. 


,33 


25.13 


29.61 


0.329 


8556 


143.38236 


+28.75308 


3. 


,33 


24.91 


29.66 


0.258 


10849 


143.36004 


+28.73414 


3. 


,35 


25.59 


29.66 


0.357 


17175 


246.46910 


+26.89244 


3. 


,34 


24.69 


29.05 


0.235 


74113 


149.88620 


+2.276064 


3. 


,33 


23.56 


29.66 


0.062 


50989 


149.83421 


+2.416729 


3, 


,31 


24.61 


29.25 


0.223 


32388 


149.77887 


+2.229502 


3. 


,30 


24.63 


29.60 


0.163 


51227 


149.89208 


+2.414816 


3. 


,28 


24.45 


29.28 


0.202 


1723 


150.44702 


+2.347633 


3. 


,30 


23.63 


28.23 


0.234 


13903 


150.41495 


+2.158999 


3. 


,29 


22.59 


28.79 


0.056 


stack 






3.3 


24.85 


30.73 


0.050 



